Inelastic scattering in a local polaron model with quadratic coupling to bosons 
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We calculate the inelastic scattering probabilities in the wide band limit of a local polaron model 
with quadratic coupling to bosons. The central object is a two-particle Green function which is 
calculated exactly using a purely algebraic approach. Compared with the usual linear interaction 
term a quadratic interaction term gives higher probabilities for inelastic scattering involving a large 
number of bosons. As an application we consider the problem hot electron mediated energy transfer 
at surfaces and use the delta self-consistent field extension of density functional theory to calculate 
and compare coupling parameters and probabilities for exciting different vibrational modes of CO 
adsorbed on a Cu(lOO) surface. 
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I. INTRODUCTION 

The local polaron model describes a localized elec- 
tronic state which is coupled to a boson field. The local 
state is then assumed to be hybridized with a continuum 
of delocalized states and are thus not an eigenstate of the 
electronic part of the Hamiltonian. 

One of the first applications of the model was the cou- 
pling of plasmons to core hole&i. and valence holes^ in 
metals. The boson field then represents the plasmons 
which can be excited by the introduction of a struc- 
tureless core-hole or a valence-hole which may be hy- 
bridized with metallic states. Plasmon excitation spec- 
tra are typically measured using Electron Energy Loss 
Spectroscopy (EELS) where the energy loss of highly en- 
ergetic electrons are measured after transmission through 
a metallic film. A similar application is that of deep-level 
spectroscopy^ where the boson field represents a phonon 
system which can be excited by the introduction of a 
core-hole. Hybridization of the core-hole is then intro- 
duced to capture the degeneracy of the core-hole with 
a continuum of states with no core-hole but a high en- 
ergy Auger electron present. A somewhat different line 
of application is that of certain rare earth compounds 
which are known to give rise to mixed valence states^*^. 
These states are characterized by an alternating valence 
in an otherwise periodic lattice which can result in un- 
usual thermodynamic properties. The reason is that the 
difference in valence results in a difference in ionic radii 
and the extra valence electron thus have a strong coupling 
to the phonon system. The model designed to capture 
the effect consists of a localized /-state (the extra valence 
electron) coupled to a continuum of delocalized electrons 
and a phonon field coupled to the /-state. AUthough, 
there is orders of magnitude differences between typical 
plasmon and phonon energies the physics in the models 
are very similar and only the model parameters differ. 

Finally the local polaron model has been applied to the 
problem of resonant tunneling^ in the context of elec- 
tronic transport, and the very similar problem of Hot 
Electron Femtochemistry at Surfaces- (HEFatS). The 



idea of HEFatS is that an adsorbate system on a metal 
surface can have unoccupied electronic states which ob- 
tain a broadening due to interaction with the metallic 
states. If a hot electron (an electron above the Fermi 
level) is generated in the metal it may interact with the 
unoccupied state and induce a chemical reaction on the 
surface. As an example we can think of a single molecule 
on a metal surface with one unoccupied electronic state 
well above the Fermi level. A hot electron with an energy 
that matches the unoccupied orbital has the possibility 
of tunnelling from the metal to the molecule resulting in 
a transient occupation of the orbital. If the molecule was 
initially in an equilibrium position the electron will as- 
sert a force on the internal molecular degrees of freedom 
and can excite vibrational modes of the molecule before 
it tunnels back into the conductor. The molecule may ac- 
quire enough energy in this process to undergo a chemical 
reaction or a desorption event. A clever method to pro- 
duce hot electrons is based on a Metal-Insulator-Metal 
(MIM) heterostructure as suggested by Gadzuk*. With 
an ideal MIM device it is possible to tune hot electrons 
to any desired resonance of an adsorbate system and the 
approach thereby suggests the highly attractive possibil- 
ity of performing selective chemistry at surfaces. Such 
devices have been constructed and characterized^ and 
comprise a promising candidate for advanced HEFatS ex- 
periments. 

Common to all these applications is that a linear 
bosonic coupling term has been assumed. It is by no 
means obvious that linear coupling captures the possibly 
complicated interaction of bosons and electrons, although 
it is probably often a good approximation. To examine 
the local polaron model beyond linear coupling we calcu- 
late the consequences of substituting the linear coupling 
term with a quadratic coupling term. In principle we 
should add the quadratic coupling on top of the linear, 
but this renders the model somewhat tedious to work 
with and the physics of quadratic coupling become hid- 
den in complicated expressions. In contrast, having only 
quadratic coupling allow us to obtain inelastic scattering 
amplitudes very similar to those with linear coupling and 
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the comparison is very instructive. In terms of bosonic 
potentials, a linear coupling term corresponds to a shift 
in the potential minimum whereas a quadratic coupling 
term corresponds to a shift in the frequency of the po- 
tential. 

The paper is organized as follows: In section [TT] we 
present the local polaron model with a general coupling 
function and no bosonic dispersion. The electronic part 
is briefly reviewed and the wide band limit which is im- 
posed in the remainder of the paper is defined. We then 
present the well known spectral function and inelastic 
scattering probabilities of the model with linear coupling 
and compare with a quadratically coupled model calcu- 
lated in the present work. It is shown that for inelastic 
scattering involving a large number of bosons, quadratic 
coupling can give rise to much larger scattering proba- 
bilities. In section Hill we apply the theory to hot elec- 
tron mediated excitation of the different modes of CO 
adsorbed on Cu(lOO). The model parameters are calcu- 
lated using density functional theory and the delta self- 
consistent field method and we find that linear coupling 
dominates desorption probabilities for the normal modes 
along the molecular axis, but vanishes for the frustrated 
rotations where quadratic coupling has to be taken into 
account. In appendix 1^ we derive a path integral rep- 
resentation of the Newns- Anderson retarded Green func- 
tion and show that the special properties of the wide band 
limit allow us to decouple bosonic and electronic degrees 
of freedom. Appendix IB] and appendix [Cl present the de- 
tails of the calculations leading to the spectral functions 
and inelastic scattering probabilities associated with lin- 
ear and quadratic coupling. In appendix |D1 we show 
how a linear transformation of creation operators makes 
it possible to obtain the exact Green functions including 
both linear and quadratic coupling. 
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II. MODEL 



The Newns- Anderson model with coupling to 
bosons 



The general model we are concerned with is composed 
of a Newns- Anderson type Hamiltonia n^^i^^ coupled to 
a dispersionless (single frequency) boson field through a 
single electronic state. A dispersionless boson field nat- 
urally corresponds to a single mode of oscillation in an 
adsorbate system, whereas we can think of the disper- 
sionless model as describing an Einstein band if the bo- 
son field represents a phonon system. Thus it is a model 
of non- interacting metallic electrons |fc), a localized reso- 
nant state I a), and a harmonic oscillator described by the 
coordinate x or equivalently the bosonic creation and an- 
nihilation operators and a. The Hamiltonian is given 



-I- hujQO^a + £a(x)c\ca, 



(1) 



where Cq creates an electron in the state \a) and Ck creates 
an electron in the state |fc). The function ea{x) couples 
the resonant electron to the oscillator degrees of freedom. 
If one considers a hole coupled to the bosons instead of 
an electron, the order of Ca and c\ should be exchanged. 

It is natural to Taylor expand the coupling function 
in the vicinity of the ground state minimum xq. Includ- 
ing only the zeroth order term £0(2^0) = results in 
the Newns- Anderson model. Since it is quadratic in the 
electronic creation and annihilation operators one could 
in principle formally diagonalize it. However, it is much 
more useful to investigate the resonant state \a) (which 
is not an eigenstate of the Hamiltonian) and we are thus 
led to consider the retarded Green function: 

Glit)^^^e{tmca{t)cl{m), (2) 

where |0) is an electronic vacuum state and 

c{t) = e'"'/^c{0)e-''"''^\ 

It is easily calculated in the energy domain using the 
Dyson equation and the result is 
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huj-eo-T.{uj)+iT{uj)/2' 



with 



and 



r(c^) = 2^^|Kfep5(;ic^-e,), 



S(^) 



duj' T{uj) 

27r LO — LO' 



(3) 



(4) 



(5) 



Assuming the hopping matrix elements Vak to be con- 
stant, r(w) becomes proportional to the metal density 
of states. If we furthermore assume the metal density 
of states to be wider than the resonance energy we can 
write r(a;) = r(£o) and E = 0. This is the wide band 
limit which will be imposed in the present paper. It allow 
us to separate electronic and bosonic degrees of freedom 
in the general case and we can calculate Green functions 
corresponding to linear and quadratic coupling exactly. 

In the wide band limit the electronic retarded Green 
function is 



G°j(t) = -i0(t)e-('^''+r/2)t/^^ 



(6) 



and the spectral function is a Lorentzian with Full Width 
at Half Maximum given by F. When boson coupling 
terms are included (first and second order Taylor expan- 
sions of £a{x) in ([T])) the spectral function changes and 
inelastic scattering on the resonance becomes possible. 
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Suppose the resonance is initially unoccupied and the 
oscillator is in the state n. The differential probability 
that an incoming particle (hole or electron) with energy 
£ will scatter through the resonance into a state of energy 
e' is given by the inelastic scattering matrix which can 
be expressed in terms of a two-particle Green function^ 
as 

i?(n;e',£) =T^ J '^^^e'ie-e-)r /r^+^e't/H-^es/r^ 

xG{n;T,s,t), (7) 

where 

G{n;T,s,t) ^ 9{s)9{t){n\ca{T - s)ci{T)ca{t)ci{0)\n). 

The probability of transferring a given amount of energy 
to the bosons can thus be obtained by integrating the 
inelastic scattering matrix over the relevant values of e'. 
In this paper we focus on inelastic scattering by elec- 
trons, since this is the relevant quantity in the context 
of HEFatS and EELS. However, the two-particle Green 
function also appear in the calculation of optical transi- 
tion amplitudes of an adsorbed moleculei^ and knowing 
G(n; T, s, t) allows one to calculate a variety of observable 
quantities. 

Finally, we note that the lifetime of the electron is 
independent of the boson coupling in the wide band limit. 
The probability that the state \a) is unoccupied and that 
the oscillator is in any state at time t given that the state 
was occupied and the oscillator was in the state \n) at 
t = is 

OO 

Pain; t)^Y.\ «; 0) I' = (8) 

m=0 

which is proved in appendix El 




■To 

X 



FIG. 1: A general example of adiabatic potentials Vi{x) and 
Vo{x) and the coupling function ea(a) = Vi{x) — ^0(2;). The 
vertical distance between the two potentials at the ground 
state minimum is eo- 




■'■0 

X 



FIG. 2: A shifted excited state corresponds to a linear cou- 
pling function Sai^) ~ x. The strength of the coupling is 
proportional to the derivative of the excited state at xq- 



B. Coupling function and adiabatic potentials 

Consider the state \x,a) with the oscillator at x and an 
electron occupying the resonance. The expectation value 
of the Hamiltonian on such a state will depend on the 
value of X due to the coupling (x) and if we could cal- 
culate the electronic energy for all values of x we would 
obtain an excited state potential Vi{x) = {x,a\H\x,a). 
Doing the same for the state with no electron in the reso- 
nance |a;, 0) would result in a different potential Vo{x) and 
the coupling function should then be given by Sa{x) — 
Vi{x) — Vo{x) which is illustrated in figure[T] In the model 
dll) we have implicitly assumed that the potential Vo{x) 
is quadratic, but in general it could have any form. The 
potentials Vi{x) and Vo{x) are called Born-Oppenheimer 
surfaces and are obtained by moving the oscillator adia- 
batically in the electronic environment. 



C. Linear coupling 

We now Taylor expand the coupling function ea{x) to 
first order and express the boson coordinate in terms of 
creation and annihilation operators. This gives an inter- 
action term: 

Hi = Xiclca{a^ +a), (9) 

with 

>^i = -^^V, , l=\[^, (10) 

V2 ox x=Xo V ITlUlQ 

where m and loq are the mass and frequency of the os- 
cillator. This model corresponds to the potentials and 
coupling function shown in figure O 
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As shown in appendix [XI the bosonic degrees of free- 
dom decouple from electronic degrees of freedom in the 
wide band limit and the retarded Green function thus 
becomes a product of an electronic part given by ^ and 
an bosonic part. Since the interacting term is linear in 
the oscillator coordinate the oscillator part of the Hamil- 
tonian can be diagonalized by "completing the square" 
Or, equivalently, performing a canonical transformation 
which shifts the boson coordinate an amount propor- 



tional to clca- 



H 



P 



Ai 



-«-^clca(a^ - a). (11) 

The retarded Green function can then be calculated ex- 
actly for the n'th excited state giving 

xL„[gi|l-e-°*n, (12) 

where L„ is the n'th Laguerre polynomial and gi = 
Xl/{huJo)'^. In this paper (;„ denotes a dimensionless effec- 
tive coupling constant and G'"-' denotes the exact Green 
function corresponding to a coupling term £a(x) ~ 
and not the contribution from an n'th order perturba- 
tive calculation as is sometimes custom. The spectral 
function is given by A'k^uj) = -2lmG^^)(n;u;) and for 
the ground state we obtain 



(13) 



m! (fiw - £o + (.91 - TO)fi^o)2 + (r/2)2 ■ 

The spectral function is thus a sum of Lorentzians of 
width r and an internal spacing of loq and the amplitude 
of the m'th peak follows a Poisson distribution. It should 
be noted that the peaks do not represent excited states of 
the oscillator. It is the spectral function of the resonant 
electron with the oscillator in the ground state and the 
different peaks show that the coupling term mixes the 
eigenstates of the isolated oscillator. The real part of 
the self energy is always negative and given by —hiuogi 
and all physical observables are invariant to Ai — Ai 
since linear coupling corresponds to a shifted harmonic 
oscillator and the direction of the shift is irrelevant. 

The two-particle Green function and inelastic scatter- 
ing matrix can also be calculated exactly^' and the 
probability that an incoming electron scatters on the res- 
onance and excites the oscillator from the ground state 
to the n'th excited state is 



pW(e) = r2e-^^4|FW(e)p, 



(14) 



with 



k=0 ^-^^ 



E 



9i 



The probability of exciting the n'th vibrational state 
thus essentially conserves the Poisson distribution, but 
the Lorentzians are replaced by the interference factor 
|F„(£)p. The results and (HI]) can also be obtained 
using a disentangling theoreroi^ as shown in appendix IB] 



D. Quadratic coupling 




FIG. 3: A frequency shifted excited state gives rise to a 
quadratic coupling function: £a{x) ~ x^. 

We will now consider an quadratic excited state po- 
tential energy surface Vi{x) which has a minimum that 
coincides with the ground state minimum, but has a dif- 
ferent harmonic evolution. The potentials and coupling 
function corresponding to this is shown in figure [3l Al- 
ternatively we could regard this model as a second order 
Taylor expansion of the phonon coupling function ea{x) 
when the first order contribution vanishes. The interac- 
tion term in the Hamiltonian becomes 



with 



A, = 



h 1 d^iVi - Vo) 



2raujQ 2 



dx"^ 



(15) 



(16) 



In this work we will only consider bound excited state 
potentials of the form Vi{x) = mujf{x — a;")^/2 and we 
can then write 



UJi = lVo^/l + Vy2JhlLJo. 



(17) 

In the wide band limit the electronic and bosonic degrees 
of freedom decouple and the boson propagator can be 
evaluated using a generalization of the Baker-Campbell- 
Hausdorff formulaF^. As a result the retarded Green 
function and spectral function can be calculated exactly. 
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The derivation is shown in appendix [C] and gives for the 
ground state: 



E 



(18) 



h 



{hLU-eo + h{ujo - wi)/2 - 2mhuji)^ + (T/2y 



where 



92 



1 9" 
ml dx^ 



-(1 



-1/2 



This result is vaUd for > which impHes that g2 < 1- 
Again, the spectral function is a sum of Lorentzians, but 
with the m'th peak damped by a factor of bm9T instead 
of a Poisson distribution. The real part of the self en- 
ergy is now given by half the frequency shift. The in- 
ternal spacing between the peaks is 2uji and we see that 
the quadratic coupling only mixes the oscillator ground 
state with the even excited states of Vi{x). This is due 
to the mirror symmetry of Sa (x) which implies that only 
oscillator states with equal parity mix. In the quadratic 
case, the effective dimensionless coupling g2 is not sim- 
ply given by X^/ifiuJo)'^ as may have been anticipated, 
and thus the m'th term in does not correspond to a 

m'th order perturbative calculation of Gj^ {t) in A2. The 
calculation leading to the exact result pE]) is very differ- 
ent from the perturbative approach and we have checked 
that the second order Taylor expansion of ^8]) indeed 
gives the result obtained from second order perturbation 
theory. 

AUthough the spectral function (|18p shows a series of 
peaks spaced by 2u!i it is only possible to excite an inte- 
ger number of cuq through inelastic scattering. The reason 
is of course that the boson field is completely decoupled 
from the electronic states in the asymptotics of a scatter- 
ing event and will thus be observed in an free oscillator 
eigenstate. Again, the symmetry of the quadratic cou- 
pling means that transitions involving an uneven number 
of bosons are forbidden. The two-particle Green function 
and inelastic scattering matrix are calculated in appendix 
[Cjand the probability for an incoming hot electron to ex- 
cite 2n quanta of oscillation when initially in the ground 
state is: 



(2) 



(19) 



with 



,.t(-.'(»)£f^^^a^ 

j=0 ^-^ ' fc=o i=0 ^ ' 

1 

e - eo + (fii^o - ftwi)/2 - 2(j + fc + l)hwi + lY 12 



The structure is very similar to the case of linear cou- 
pling. With linear coupling the probability for an elec- 
tron to create n bosons are proportional to the n'th order 
Taylor expansion of e^^ and normalized by e~^^^ whereas 




2 4 6 

(e-ej/wo 



2 4 6 

(e-ej/^o 



FIG. 4: Probabilities of making the transition — > n through 
resonant inelastic scattering with linear and quadratic cou- 
pling. The parameters are F/fttJo = 0.5, g\ = 0.2, and 
LO\ = 0.75a;o (g2 = 0.02). Even though 32 < 5i the quadratic 
coupling becomes dominating for large n due to the slowly 
decaying expansion coefficients. One should also note the 
spacing between peaks which is loq for linear coupling and 
2lji for quadratic coupling. The centers of the probability 
distributions are approximately shifted by nLOo/2 for the lin- 
ear coupiing and naji/2 for the quadratic coupling relative to 
the bare resonance energy eo , since this is where the binomial 
coefficients in (|14|l and (|19p have their maxima. 



in the quadratic case the probability to create 2n bosons 
are proportional to the n'th order Taylor expansion of 
(1 — 52)"^^^ and normalized by (1 — 52)- In the context 
of EELS and plasmon excitations, one would now ob- 
serve a series of peaks spaced by 2u}q. If the plasma fre- 
quency is not known the spacing itself cannot give clues 
to whether linear or quadratic coupling governs the tran- 
sitions, but one could use the relative amplitude between 
peaks since these follow a Poisson distribution if linear 
coupling dominates and the distribution bng2 if quadratic 
coupling dominates. If both linear and quadratic cou- 
pling is present one would observe a coupling dependent 
enhancement of every second peak. 

In a model with linear coupling, the probability of 
exciting 2n vibrational quanta is proportional to 5^" 
whereas it is proportional to in a quadratic coupled 
model. This implies that if (72 > 5i a quadratic coupling 
term will give rise to larger inelastic scattering probabil- 
ities than a linear term. Even with 52 < 5i a quadratic 
coupling term may have stronger effect for large n since 
the expansion coefficients of {l — x)~^/'^ decay slower than 
those of e^. This is illustrated in figurelH where the prob- 
ability of transferring n vibrational quanta to the ground 
state is shown for linear and quadratic coupling. 
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III. APPLICATION TO HOT ELECTRON 
MEDIATED DESORPTION 

As an example of a system where the dynamics can be 
approximated by a local polaron model, we consider the 
problem of hot electron mediated energy transfer on a 
metal surface. Such an energy transfer can lead to des- 
orption of adsorbed molecule9i^>^'^»^'^>^i2ii22, or in- 
duce chemical reactions which cannot proceed by thermal 
heating. The conceptual picture of the process is the 
following: Hot electrons are generated in the metal by 
means of an MIM device or a femtosecond laser. The hot 
electrons may then interact with a chemisorped molecule 
by tunneling from the metal to an unoccupied molecu- 
lar state and excite vibrational states in the molecule. If 
enough energy is transferred to the molecule either by 
a single or multiple scattering events, the molecule may 
eventually desorp or break an internal chemical bond. As 
a particular example we will calculate transition proba- 
bilities for CO adsorbed on Cu(lOO). 

To calculate inelastic scattering probabilities within 
the local polaron model we need to obtain the coupling 
fmiction Sa {x) . As described in section III Bl we can fix 
the molecule at different positions and calculate the po- 
tential energy surfaces Vi{x) and Vq{x) at each point 
and ea{x) = Vi{x) — Vo{x). The model ([1]) does not 
directly contain Coulomb interactions between electrons, 
but these are included in the calculation of Sa {x) which 
thus becomes an effective coupling that is supposed to 
contain all the electronic interaction associated with the 
excited state of the molecule. 

The potential energies Vi{x) and Vb(x) has been ob- 
tained using the code gpawS^i^ which is a real-space Den- 
sity Functional Theory (DFT) code that uses the projec- 
tor augmented wave 

method.26'27 aU our calculations 
we used the Revised Perdew-Burke-Ernzerhof (RPBE) 
exchange-correlation functional2^ since this has been de- 
signed to perform well for molecules adsorbed on surfaces, 
and has been shown to perform better than the original 
PBE functional^^ for adsorbed molecules. 

We set up a Cu(lOO) surface consisting of three atomic 
layers with the top layer being relaxed. 10 A of vac- 
uum has then been introduced above the slab and 0.50 
monolayer of adsorbate molecules relaxed at top sites 
which is the preferred adsorption site. Both molecules 
adsorb with their molecular axis perpendicular to the 
surface with O pointing away from the surface. We 
then did a normal mode analysis and mapped out the 
three ground state potential energy functions Vo(a;i) cor- 
responding to the two normal modes that involve the 
perpendicular degrees of freedom and a frustrated rota- 
tion. The perpendicular modes roughly correspond to 
an internal stretch d = xq — Xc(n) s-nd center of mass 
z = {moxo + mc(N)Xc(N))/ {mo + mciN))- We do not 
include the three remaining molecular modes since one is 
another frustrated rotation with identical properties to 
the one considered, and the two frustrated translations 
are only weakly coupled to the resonant electron and are 



CO on Cu(100) 



^ -69.0 




-0.4 -0.2 0.0 0,2 0,4 

X - 1(1 [A] 



FIG. 5: Potential energy surfaces along the frustrated rotation 
mode of CO adsorbed on a Cu(lOO) surface. The coordinate 
a; is a generalized coordinate representing the deviation from 
equilibrium, x = 0.4 corresponds to a 24° angular deviation 
from the perpendicular position. 



not expected to play a significant role in the femtochem- 
istry. In all calculations we use a p(2x2) cell, sample 12 
irreducible k-points in the surface plane, and use a grid 
spacing of 0.2 A. 

To find the excited state potential energies Vi{xi) cor- 
responding to the three normal modes of interest, we have 
used the method of linear expansion ASCF which has 
been published in a previous work'^S and implemented 
in gpaw. In the previous publication we have tested the 
method against inverse photo-emission spectroscopy, and 
found that it performed well for molecules chemisorped 
on surfaces In each step of the self consistency cycle 
an electron is removed from the Fermi level, the density 
of an excited state is added to the total density, and the 
band energy of this state is added to the total energy. 
To get the band energy right we need to expand the ex- 
cited state on the Kohn-Sham orbitals found in each it- 
eration. The method is thus a generalization of the usual 
ASCF where occupation numbers are changed. Instead 
of changing occupation numbers we occupy an orbital 
which is not an eigenstate of the Kohn-Sham Hamilto- 
nian but a superposition of eigenstates in such a way that 
the state is as close as possible to the original molecu- 
lar state. In the present case the excited state is the 
anti-bonding 27r orbital of CO. In figure [5] we show the 
ground and excited state potential energy surfaces corre- 
sponding to the frustrated rotation along with Saix). It 
is clear that the excited state potential is not exactly a 
quadratic potential and the parameter A2 which we need 
to calculate transition probabilities will depend on how 
we fit this potential to a quadratic form. However the 
width of the Gaussian ground state vibrational wave- 
function corresponds to x = 0.08 A and for low lying 
excitations we can thus use this region of the potential 
which is rather flat. In fact, a closer look at the excited 
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Mode 


hu) 


Ai 


Aa 


Frustrated rotation 


0.037 





-0.009 


Center of mass 


0.043 


-0.006 


~ 


Internal stretch 


0.248 


-0.170 


~ 



TABLE I: Parameters for CO adsorbed on Cu(lOO). All num- 
ber are eV. Note that while the quadratic coupling for the 
two perpendicular modes are very small and thus neglectable, 
the linear coupling of frustrated rotation vanishes exactly due 
to symmetry. 



state potential reveals that the ground state minimum 
geometry actually has an unstable extremum in the ex- 
cited state, but since the curvature is rather small we 
will simply approximate it by a constant potential. For 
both perpendicular modes we find that ea{xi) ~ Xi and 
quadratic coupling can thus be neglected. In contrast, 
due to symmetry the excited state potential energy of 
frustrated rotation is invariant to Xi —Xi and the lin- 
ear coupling term thus vanishes. We have calculated the 
excitation energy to Saixa) = 2.8 eV and the resonance 
width is estimated from the Kohn-Sham projected den- 
sity of states to r « 1.0 eV. In table HI we display the 
calculated parameters corresponding to the three modes. 

We note that when calculating transition probabilities 
we should include all modes in the model ([1]), because 
even if the modes are not coupled directly they have an 
indirect coupling since they all interact with the reso- 
nance. It is possible to obtain expressions for the scat- 
tering matrix including more than one mode, but these 
are rather complicated to handle and for weakly cou- 
pled systems the physics can usually be extracted from 
three one-mode models.^^^ In figure [6] we show the calcu- 
lated probabilities for a hot electron to excite the different 
modes of CO adsorbed on Cu(lOO). The internal stretch 
and and frustrated rotation show transition probabilities 
on the same order of magnitude whereas the center of 
mass vibrations are very unlikely to get excited. This is 
in accord with calculations of the electronic friction co- 
efficients of this systero^i^^ which is very closely related 
to the coupling function (x) M The frequency of inter- 
nal vibration is five times larger than both the center 
of mass and frustrated rotation frequencies and as previ- 
ously showrti^ the stretch mode will completely dominate 
the total energy transfer. Thus, in a simple model where 
hot electron mediated desorption^^ is reduced to calcu- 
lating the probability of transferring the chemisorption 
energy to the adsorbate, the internal mode governs the 
desorption probability. Nevertheless, our estimate of F 
is based on the Kohn-Sham density of states which may 
give a poor description of the electronic spectral function 
A^{u!). If F is significantly smaller than our estimate, the 
quadratically coupled frustrated rotation will play an im- 
portant role in hot electron mediated desorption for this 
system. 



IV. SUMMARY AND DISCUSSION 

We have calculated the spectral function and inelas- 
tic scattering amplitudes in a local polaron model with 
quadratic coupling to bosons. The probability of exciting 
n bosons is found to be damped by a distribution func- 
tion given by the n'th Taylor expansion of 1/^/1 — g2 
which decays much slower than the Poisson distribution 
appearing in a linearly coupled model. Hence for compa- 
rable values of linear and quadratic coupling constants, 
a quadratic term will dominate inelastic scattering prob- 
abilities involving a large number of bosonic excitations. 

As an application we have considered the prob- 
lem of hot electron mediated vibrational excitations of 
molecules adsorbed on metal surfaces. The coupling con- 
stants were calculated from the excitation energy along 
the molecular normal modes using delta self-consistent 
field DFT. It was found that quadratic coupling is im- 
portant for exciting the frustrated rotations since this 
mode does not couple linearly due to symmetry. 

A major approximation in the model is the quadratic 
assumption for the ground state potential. In our numer- 
ical example with HEFatS it is clear from figure [5] that 
the potentials is not exactly quadratic. For the center of 
mass mode the situation is even worse and a Morse po- 
tential is much better suited to describe this mode. The 
anharmonic deviations are likely to have a significant ef- 
fect on high lying excited states, but renders the model 
much more complicated. In fact, since the coupling to 
the internal stretch mode seems to govern the rate of 
energy transfer, one has to assume that the energy is 
readily redistributed to other degrees of freedom and an- 
harmonic coupling is thus expected to play a vital role in 
the actual desorption process. 

The wide band limit has been essential in the deriva- 




FIG. 6: Probabilities of exciting two and four quanta of vi- 
brations to the center of mass, internal stretch and frustrated 
rotation modes CO adsorbed on Cu(lOO). 
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tion of scattering amplitudes, and we do not have the 
means to solve the model (with linear or quadratic cou- 
pling) exactly beyond this approximation. However, it 
would be very interesting to do perturbation theory with 
the general retarded Green function ([3]) to examine the 
effect of energy dependence in the electronic self energy. 
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APPENDIX A: DECOUPLING OF ELECTRONIC AND BOSONIC DEGREES OF FREEDOM IN THE 

WIDE BAND LIMIT 



1. Path integral representation of the Newns- Anderson retarded Green function 



The path integral representation of propagators often give a renewed insight into the underlying physics, allthough 
the mathematical complexity can be somewhat larger. By writing the Newns-Anderson retarded Green function as 
a sum over paths, we see that each path can be understood as sequence of jumps from the resonance to the metallic 
band and we then have to sum over all possible time intervals between each jump. In the wide band limit the time 
spend in the metal band goes to zero and the electron thus spends all the time of propagation in the resonant state. 
The Newns-Anderson model is given by ([!]) with ea{x) = Sa{xQ) — Eq: 

Ho = eocica + J2 ^kcick + J2 {Vakcick + Kfe4ca) (Al) 

k k 



and 



G%it) = -^eit){a\e'^"^a). 



(A2) 



The path integral representation is derived by dividing the time interval t in N intervals of length At — t/N. When 
N becomes sufficiently large we can take e~*^* = (^g-'^HAt-jN ^ _ m/^t^N ^ -yy^ then insert iV — 1 complete sets of 
states \n) such that the Green function becomes a sum over iV-fold products of matrix elements 

G°R{t) f=i -iO{t) ^ {a\l-iHQM\ni){ni\l-iHM\n2) ...{nN^i\l-iHo/^t\a). (A3) 

ni,n2---njv-l 



Assuming that (a|fc) = the states \n) can either be \a) or \k) and the matrix elements (a|l — iHo/S.t\a) — e^**^""^*, 
(/cijl — ii/oAt|fc2) — (Jfeife^e""'"^* and {a\\ — iHo/S.t\k) = —iVak^t represent propagation in the resonance, propagation 
in the band and a jump from band to the resonance respectively. When we take the limit TV — > oo, (jASp becomes 
formally exact and the jumps between band and resonance become instantaneous. It is then most convenient to order 
the terms in (jASp according to the number of jumps. Since the endpoints of the time interval is at the resonance, a 
jump into the band has to be accompanied by a jump back into the resonance and each such "band excursion" comes 
with a factor of — l^afepe~*''='^' where is the time spend in the z'th excursion into band. It is also clear that p 
excursions into the resonance has to be accompanied by p + I resonant propagation factors e~''^°''^ where ai is the 
I'th time interval in the resonant state. Finally, for a given number of band excursion we have to integrate over all 
possible band and resonance time intervals and the retarded Green function becomes 



poo / poc />oo \ 1 



X(5(ao + 5I(^j+^j)-0' (A4) 
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where the deha function has been introduced to ensure that the time intervals sum to t and the theta function has 
become redundant. We can use the delta function to eliminate the ai integration variables and get 

//•OO / pOG 1*00 \P 

p^e^5o(x)=_,e-^eoty^ da^Y^y- dadTT{T)y(ao + Y.{a,+T,)-t), (A5) 



with 



/OO T 
-^^^r(c^)e-^("-^°)*. (A6) 
-00 27r 

By using that 5{t) = (l/27r) J dwe*'^* it is now possible to evaluate (|A5p and recover the result ([3]). In the wide 
band limit, T[t) — T5{t), which implies that the electron does not spend any time in the band and the retarded 
Green function becomes a sum over paths which are composed of instantaneous excursions into the band. We use the 
notation x to represent a position in state space and J T>xe'^^° ^"^^ as a formal expression representing the sum over all 
paths weighted by the Newns- Anderson action Sq. 



2. Resonant electron in a bosonic environment 



We now proceed with the full Hamiltonian (jT]). Introducing the bosonic coordinate x the full retarded Green 
bmction (with the boson field in the state a;o) can be written 

Gr (xo ; = y VxVxe'^^ (^)+*^« ('^'^^ , ( A7) 

where 5*0 (x) is the Newns- Anderson action, Sb{x) is the free bosonic action corresponding to the Hamiltonian Hb = 
ujQa'^a, and Sj is the interaction part of the action corresponding to the Hamiltonian Hj = c^CaSaix). However, 
there is a much nicer way to handle the coupling to the boson field. One can think of the bosons as an environment 
influencing the paths of the resonant electron and it can be shown that the Green function can be written^i^^ 

GR{xo;t) = J Vxe''"^^Hxo\U{x;t)\xo), (A8) 

with the environment time evolution operator 

U{x;t) = e^^'^'Te-'^odt'Hiixit')]^ (A9) 

where Hj{x{t')) — Hb + Hj{x{t')) is the environment Hamiltonian evaluated on an electronic state fixed at x and T 
denotes time ordering. Thus, the price we pay in separating bosonic and electronic degrees of freedom is an explicit 
path dependence in the environment part of the propagator. In general it is not possible to evaluate Hj{x{t')) on all 
possible paths, but in the wide band limit it is particularly simple. The reason is that the resonant electron stays on 
the resonance in all possible paths and the environment Hamiltonian is therefore independent of the electronic path. 
In fact, the environment propagator becomes U{x',t) = g*'^a+atg-jwotot-jea(a:)t ^^iq electronic and bosonic degrees 
of freedom completely decouple in the retarded Green function: 

GRin;t) = Glit)GBin;t), GB{n;t) = („|e*'^«*«*e-'"'^^'^*-^'^-(^)*|n). (AlO) 

The situation is very similar for the two-particle Green function: 

G{n; r, s, t) = e{s)e{t){n\ca{T ~ s)4(r)c4i)4(0)|n) = e{s)e{t){n\ca{T - s)U{t - s, T)ci{T)ca{t)U{t, 0)ct (0)|n). 

The resonant electron is first propagated forward in time from to t and then backward in time from t to r — s. 
The interaction vanishes between the c^(t) and c{t) because the resonant state is unoccupied here. Again, the full 
Green function can be written in terms of a bosonic influence propagator and in the wide band limit the electronic 
and bosonic degrees of freedom decouple so: 

G{n- T, s, t) = Gl{t)Gl{s)GB{n; r, s, t) (All) 
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with 

Gb(71; T, S, t) = {nyHo{r-s)^^(H^+e^(x))s^-^H^r^^H^t^-^(H^+e^(x))t^^y (^^2) 

This can be seen by applying the arguments above to both of the time evolution operators and the fact that G°(t, s,t) — 
G°«(t)G?j(5). 

We observe that for any coupling function Sa{x) the above form of the two-particle Green function implies that 
G{n;t,t,t) = |G5j.(t)p. This means that in the wide band limit, the resonant lifetime is unaffected by the phonon 
coupling since the probability of finding the electron in the state \a) at time t (and the oscillator in any state |m)) 
given that it was there at i = (where the oscillator was in the state |n)) is 



Pa{n;t) = ^ \{m,a;t\n,a;0)\'^ = ^ \{m\ca{t)cl\n)\'^ = {n\cacl{t)ca{t)cl\n) 

m— m— 

= G{n- 1, t, t) = |G?j(t)p = e^". (A13) 



So, in the wide band limit the resonant state always has a well defined lifetime given by Ta — h/T 

The problem of calculating the inelastic scattering matrix has now been reduced to evaluating the phonon propagator 
Gsin; T, s, t). In general this is not an easy task, but we will show that in the case of linear and quadratic coupling 
terms we can use a disentangling theoremi^ to write the exponential operators in a form that allows a direct evaluation 
of the expectation value. The theorem is a generalization of the Baker-Campbell-Hausdorff theorem and states that 
if A, B, and G are three operators with a closed commutator algebra then QO-A+bB+cC _ j^^aA^pB^jC ^ where h, a, /?, 
and 7 are known functions of a, b, c, and the commutation parameters. 

APPENDIX B: GREEN FUNCTIONS WITH LINEAR COUPLING TO BOSONS 

1. Single-particle Green function 

The model ([T]) with linear coupling function given by 

ea{x)^Xi{a^ +a), (Bl) 

is well-known and the one-particle Green functions can be obtained exactly in the wide band limit by a canonical 
transformation.i^^ Here we derive it using the formalism above and the disentangling theorem.— To obtain the 
one-particle Green function we need to evaluate the boson propagator 

G^B\n-:t) = (^|e^"o*°*°e-''^°*(''^''+^i(°+"^^/""^|7i). (B2) 
Using the disentangling theorem on the second exponential operator leads directly to the expression 

where Ln{x) is the n'th Laguerre polynomial. To obtain the ground state spectral function p3p we Taylor expand 
exp(gie"^''*), perform a Fourier transformation, and take the imaginary part. 

2. Two-particle Green function 

The procedure can also be used to obtain the two-particle Green function. The object of interest is now the 
two-particle boson propagator which we write 



q{^) Y s t) — ^^|g*'^o('r-s)a^"agMos(ata+Ai(a+at)/cJo)g-i'^oT'a*agii^o*a^ag-'ii^ot(a*a+Ai(a+at)/cJo) 
_^inwo{'r-t)^-gi{l-iiAJot-e^'-'^o') ^-gi{l+iuJos-e"^o") 

x(n|e'^o^ ' e '^o^ ' g ^-"^ "g -^o^ ' g-^o^ ' (B4) 
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where we used two different forms of the disentangUng theorem to move e"^"*"^" and e ^'^ota^a ^j^g ^eft and right 
respectively. We then use the theorem again to move e^^ft*""^)"^" to the left and obtain 

Gg^(„; r, S, t) ^e-9i(2-»-o(t-.)-e-'^o*-e-o.) 

FinaUy, we can use the Baker-Campbell-Hausdorff theorem to collect all lowering operators at the the right. The 
evaluation of the remainder gives a Laguerre polynomial and the result is 

G%\n-T,s,t) = e's^-«(*-^)e-«^/-o(-^^.*)L„([g,(/^^ (B5) 

with 

/^^(r, s, i) = 2 - e-"^"' - e"^°' + e-"^°''{l - e*'^°*)(l - e"^°'). (B6) 
The inelastic scattering matrix has previously been derived^ii^ and we will not repeat the calculation here. 



APPENDIX C: GREEN FUNCTIONS WITH QUADRATIC COUPLING TO BOSONS 

1. Single- particle Green function 

We now consider a quadratic coupling function of the form 

ea{x) ^ X2{a^ + a)^ (CI) 
To obtain the retarded Green function we would like to calculate the boson propagator 

G^'^\n' t) = ^^^^i'^ota'' a ^-iwoa'^ at-iX2ia^ a+aa'* +aa+a^ )t^^\^ _ gniwo* g-'jA2 1 ^^^j g-i(wo+2 A2)at at-iA2 (aa+a+a+)t (C2) 

We proceed by disentangling the exponential operator: 

g-i(i^0+2A2)atat-A2(aa+a+ot)t _ ^i\2t ^iujo /2t /2 ^fa^ ^fe^^" aa ^ga^ a (C3) 

where 

-A2tanh(iwit) / + 2A2 . ,. \ , , 

fit) = — — g{t) - In cosh(iwii) + smh(iwit , (C4) 

LOi + (cjo + 2A2j tanh(^a;l^j V wi / 

and LOi — Li;o(l + 4A2/clIo)^/^. This is valid for a bound excited state potential with positive second derivative in which 
case the argument of the square root is positive. In the case of an unbound excited state potential, the functions / 
and g involve real hyperbolic functions and the spectral function acquires a qualitatively different structure. Acting 
with the operator e""" on a state \n) gives 

[»/2] I , 1 

^""»=E|((;^)^l-20, (C5) 

where [n/2] means the integer part of n/2. Collecting it all and noting that /(i)e~^^*^ — sinh(ici;it), gives the 
retarded Green function in the wide band limit: 

OA 1 [^/2] J 2/ I 

Gf{n,t) = -^0(^)e(-«-^/^)*e'("+l/2)^ot^osh(^c.l^) + ^^^i±^ sinh(^c.it) 7Ju2r^^ 

V LOi / [I.) (n 21). 

with 

A2 

h= — sinh{iujit) (C6) 
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To find the spectral function of tlie oscillator ground state we Taylor expand the square root and obtain 



0^0 — t^i + 2A2 _ 
Wo + ^^1 + 2A2 



2wi 



1^ 00 

m=0 



with 



52 



CiJo — ti^l + 2A2 (LlIq — UJ\ 



CJQ + t^l + 2A2 

Fourier transforming and taking the imaginary part then gives 

Omi/2 



1 a™ 

m! dx" 



-(\-x) 



-1/2 



cji > (C7) 



(C8) 



(t^ - eo + (t^o - ^i)/2 - 2mwi)2 + (r/2)2 ' 



wi > 



(C9) 



where we also used that 2a;i/(wi + wo + 2A2) = 1 — 32- Note that the condition of wi > implies that 52 < 1- 



2. Two-particle Green function 



We now need the propagator 

^(2)^^. Y s t) — ^f^|g*'^oa*a('r — s)gi(cJo+2A2)atas+iA2(aa+ata'f)Sg-icJoa*aTgMoa'^atg— i([»Jo+2A2)a*at— iA2(aa+a'^a*)t|^^ 

oc 

— ginuJo{T-s) ^imuJa(t-T) ^^|gi(wo+2A2)atas+iA2 (aa+atat)s |g-i(wo+2A2)at at-iA2 (aa+a+at )t (CIO) 



m=0 



We restrict the calculation to the ground state two-particle Green function which using the disentangled expression 
(IC3p becomes 

m=0 ^ 

It is also possible to obtain a closed expression that does not involve the infinite sum, since instead of inserting a 
complete set in (jC10|) we could have brought all lowering operators to the left by repeated use of the disentangling 
theorem. However, to calculate the inelastic scattering matrix ([7]) we need to integrate over r which is more convenient 
in the present form. Performing the r integral and using the resulting delta function to replace 2mujo with e — e', 
leaves the two remaining integrals as complex conjugates. We note that bm = (2m)!/4™(m!)2 and write 

00 

i?(2)(£',e) = T'Y1 4'"6m'5(£ - £' - 2mu;o)\D^ie)\^ (C12) 

771—1 

with 



'dte-(-— /2-^/^)*f cosh(.c.it) + ^^^i±^sinh(zwri))"'^Y tanh(za;it) 

V ^ ^ tJi ^ V V^i + (cjo + 2A2)tanh(icJit) 



The reason we have excluded the m — term is that it does not give rise to inelastic scattering and the elastic part of 
the scattering matrix have an additional term that we do not consider here.^ This expression implies that quadratic 
coupling can only give rise to inelastic scattering events involving an even number of vibrational quanta. This is also 
true if the initial state is not the ground state, since from (jClOp we see that in general (m — n) has to be even. To 
evaluate the inelastic scattering matrix we note that 

/ -A2tanh(^Wlt) A" ^ / -A2 A"/^ p-^^^^'V ( ^ V 

V^i + (wo + 2A2) tanh(iLJii)/ V cji + cjq + 2A2 / V / V 1 - ^26"^'"^* ^ 



m/2 m 

92 — 



j=0 ^ ^ fe=0 ^ ^ 
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where we used that —\2/{loi + cjq + 2A2) = J^I2. The Taylor expansion of the square root gives 



cosh(iL^ii) + !^^^±^ sinh(iL^it) = yr^y 6;gf,e-*(2'+i/^)"^*, c^i > 0. 



leading to 



. /I m/2 m / s 00 00 



6ig^+'(m + /c- 1)! 



, 7 / — ' — ' A:!(m — 1)! 



e - eo + (wo - t^i)/2 - 2(j + fc + O^^i + ir/2 ' 



and 



i?(2)(e',e) = r2(l - .92) y 6„,g2"<5(£ - e' - 2mc^o) 



(C13) 



(C14) 



%2 + 



1 



APPENDIX D: LINEAR AND QUADRATIC COUPLING COMBINED 



It is in principle straightforward to generalize the expressions above to the case of a linear and a quadratic coupling 
term in the Haniiltonian. The linear term can be transformed away by noting that 



+ Ai(a^ + a) + A2(a^ + a)^ = ujQO^a + A2(a^ + of' ~ 7A1, 



(Dl) 



with 



a = a + 7, = + 7, 
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Ai 



Wo + 4A2 



a^l = 1. 



(D2) 



Since the commutator algebra of a and is identical to that of a and we can immediately write down the 
one-particle boson propagator in its disentangled form 



(D3) 



Using that 



we can evaluate the propagator in the ground state and obtain 



(D4) 



(D5) 



The expression clearly reduces to G^^* (t) and G^^* (t) in the limits Ai ^ and A2 — > respectively. It should be 
straightforward to obtain the spectral function by Fourier transforming this expression after a Taylor expansion of 
the exponentials. However, the result becomes rather involved and we will not attempt to do the calculation here. 
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